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ABSTRACT 

This  thesis  undertook  a  macroscopic,  theoretical  analysis  of  a  long, 
circular,  solid  propellant,  rocket  motor  undergoing  thermal  loading.   It 
was  shown  that  classical  plasticity's  von  Mises  yield  condition  adequately 
described  dewetting  of  solid  propellant  which  was  modeled  to  have  linear 
strain  hardening  characteristics.   The  mechanical  properties  of  the 
propellant  permitted  a  closed-form  solution  for  the  stresses  and  the 
strains  resulting  from  cooling  the  motor  from  160°F  to  -60°F. 
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I.   INTRODUCTION 

The  problem  analyzed  in  this  thesis  was  that  of  a  long,  cylindrical, 
solid  propellant  rocket  motor  (see  Figure  1)  containing  a  circular  bore 
and  with  a  thin  steel  case  bonded  to  its  outer  radius.   The  motor  was 
thermally  loaded  by  cooling  it  from  160°F  to  -60°F,  which  represented 
the  extremes  that  might  realistically  be  expected  during  its  lifetime, 
from  cure  at  the  upper  temperature  to  storage  at  the  lower  one.   The 
cool-down  rate  was  assumed  to  be  slow  enough  to  maintain  equilibrium 
conditions  throughout  the  propellant;  i.  e.,  temperature  did  not  vary  with 
radius.   The  stresses  in  the  motor  resulted  from  the  difference  in  the 
coefficients  of  thermal  expansion  between  the  case  and  the  propellant. 
Because  the  propellant  had  the  higher  of  the  two  coefficients,  it  tended 
to  contract  more  than  the  case;  thus  an  applied  load  was  induced  at  the 
outer  boundary  of  the  propellant. 

The  propellant  was  considered  to  be  an  isotropic,  solid  material; 
although  in  reality  it  was  a  composite  of  solid  oxidizer  particles 
imbedded  in  and  bonded  to  a  carbon  fuel  matrix.   After  an  initial  linear 
region,  nonlinear  stress-strain  behavior  under  loading  resulted  from  the 
progressive  failure  of  these  bonds  at  the  interfaces  between  the  fuel 
matrix  and  the  oxidizer  particles.   Voids  or  vacuoles  were  thus  created 
throughout  the  material,  and  their  volume  increased  with  further  loading. 
This  phenomenon  is  called  dewetting. 

Dewetting  causes  nonlinear  material  behavior  and  has  many  similari- 
ties with  yielding  in  classical  plasticity  theory,  as  was  pointed  out 


by  Lindsey  and  co-workers,  (4)  who  developed  an  isotropic  theory  for 

dewettable  solids.   In  it,  the  strain  increment  was  decomposed  into  an 
elastic  component  and  a  dewetted  component. 

d£4j  -  dE-j  +  d£;j  (1) 

They  further  developed  from  Drucker's  postulate  that  the  dewetted  strain 
increment  was  normal  to  the  dewetting  surface. 

d£:     x  ^§fe.«0  (2) 

where  f=o  describes  the  dewetting  surface.   From  this  theory,  Wood  (7) 


^h-.it^v*0,..  o) 


showed  that  the  general  form  of  f  was 

-  B   for  I,  «( 
+AI*  -B  for  3,  >  0 

where  J2  was  the  second  deviatoric  stress  invariant,  Ij_  was  the  first 
stress  invariant,  and  A  and  B  were  mechanical  properties.   Wood  experi- 
mentally measured  these  parameters  for  a  PBAN  propellant  and  found  that 
A  =  -0.0026  and  B  =  180. 

Lindsey  and  Wood  (A)  described  the  dewetting  surface  as  a  circular 
cylinder  about  the  hydrostatic  axis  in  principle  stress  space.   The 

positive  end  of  this  cylinder  was  capped  by  a  body  of  revolution  (see 

2 
Figure  2) .   Using  a  space  defined  by  an  axis  system  consisting  of  I  , 

which  was  coincident  with  the  hydrostatic  axis  and  J  which  was  perpen- 
dicular to  the  hydrostatic  axis,  Wood's  dewetting  criterion  became 
intersecting  straight  lines  with  an  extremely  small  negative  slope  of 
-0.0026  and  an  intersection  with  the  J«  axis  of  180  ,  as  shown  in 
Figure  3.   Taking  the  elastic  propellant  properties  of  Wood's  material  and 
the  elastic  solution  for  an  encased  cylinder  subjected  to  thermal  cool- 
down,  the  intersection  between  the  thermally  induced  stress  field  and 
the  dewetting  surface  was  determined  by  calculating  the  ratio  of  J2  and 
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and  I   ,  and  finding  its  intersection  with  the  dewetting  surface. 

(See  Figure  1)   It  was  found  that  only  1  1/3%  error  resulted  in  assuming 
the  dewetted  surface  was  a  cylinder  for  the  thermal  problem.   This  was 
quite  acceptable,  especially  in  view  of  the  great  mathematical  simpli- 
fication which  resulted.   Thus,  Wood's  dewetting  criterion  for  this 
problem  was  approximated  to  be 

J2  %  B(T)  (A) 

where  for  the  thermal  problem,  B  is  a  function  of  temperature.   It  was 
observed  that  this  dewetting  criterion  was  identical  to  the  von  Mises 
yield  criterion  of  classical  plasticity. 

The  goal  of  this  thesis,  then  was  to  perform  a  stress  and  strain 
analysis  using  Wood's  dewetting  criterion  as  simplified  accounting  for 
the  nonlinear  behavior  caused  by  dewetting,  on  a  long,  circular,  solid 
propellant  rocket  motor. 

Prior  theory  had  considered  propellant  to  be  completely  elastic; 
however,  the  phenomenon  of  dewetting  was  expected  .to  occur  at  the 
inner  bore  initially  and  to  expand  outward  as  loading  was  increased. 
After  yielding  occurs  in  either  dewetting  or  plasticity  theory,  the 
linear  stress-strain  law  is  no  longer  valid  and  must  be  replaced  by  a 
flow  rule.   Both  theories  insist  upon  a  similar  form  of  the  flow  rule 
and  for  the  dewetting  surface  of  equation  (4) ,  the  flow  law  of  equation 
(2)  becomes 

where  d£;.  is  the  increment  of  dewetted  strain,  S.--  is  the  deviatoric 
stress  component,  and  a  is  a  scalar  incorporating  the  loading  history. 
Mendelson  (5)  shows  how  K  can  be  written  in  terms  of  otf  an  equivalent 
stress,  and  o|£pan  equivalent  dewetting  strain  increment: 

A    -  l4h  (6) 
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In  this  analysis  both  the  dewetting  function  and  the  propellant's 
mechanical  properties  (except  for  Poisson  ratio  which  was  assumed 
constant)  were  functions  of  temperature.   These  properties  were  determined 
experimentally  for  a  PBAN  propellant  manufactured  by  United  Technology 
Center,  Sunnyvale,  California.   Thus  in  the  solution  as  the  temperature 
loading  was  changed,  it  was  possible  to  account  for  variable  dewetting 
requirements  as  well  as  variable  stress-strain  relationships. 


II.   BACKGROUND 

Because  of  the  close  analogy  between  dewetting  and  plasticity,  the 
pertinent  literature  of  plasticity  was  searched  for  solutions  to  cylinder 
problems.   It  was  found  that  the  subject  of  partially  plastic,  thick- 
walled  tubes  received  much  attention  from  many  authors  during  the  late 
forties  and  the  fifties.   In  general,  the  vast  majority  of  the  work 
dealt  with  internally  pressurized  cylinders,  which  were  then  of  current 
interest  and  were  more  amenable  to  solution  than  those  problems  which 
deal  with  thermal  loading. 

In  1957,  Steele  (6)  made  a  review  of  several  theories  for  partially 

plastic  thick-walled  tubes.   He  stated: 

"The  mathematician's  approach  to  the  problem  under 
consideration  is  fundamentally  different  from  that  of 
the  engineer.   The  former  is  interested  in  an  accurate 
knowledge  of  the  stresses  and  strains  in  terms  of  the 
applied  load  or  loads  for  all  points  in  an  idealized 
body.   In  general,  the  engineer  seeks  an  estimation  of 
the  variation  of  factors  with  applied  loads,  which 
are  of  significance  only  in  the  functional  design  of 
the  part....  In  the  case  of  the  thick-walled  cylinder, 
interest  lies  in  the  deflections  of  the  bore  and  the 
outer  surface  for  specific  internal  pressures." 

As  an  example,  Hill,  Lee  and  Tupper  (2),  in  a  rigorous  work,  showed 
that  Nadai's  approximate  solution  for  the  axial  stress  was  in  error  by 
sixty  percent,  although  the  deflections  at  the  bore  and  outer  surface 
correlated  closely.   Rigorous  solutions  such  as  Hill's,  Lee's  and 
Tupper 's,  did  not  yield  a  closed  form  solution;  and,  although  mathema- 
tically correct,  they  were  so  complicated  that  their  usefulness  to  the 
engineer  was  minimal. 


A  general  outline  of  the  relationships  required  by  the  thick-walled 
cylinder  problem  is  given  below.   These  are  required  for  any  mode  of 
loading,  and  apply  to  elastic,  plastic,  or  dewetting  theories  equally. 

(1)  Radial  Equilibrium 

(2)  Strain  Compatability 

(3)  Stress-strain  Relation 

(4)  Continuity  of  radial  stresses  and  displacements  across  the  yield 
surface 

(5)  Boundary  Conditions 

(6)  Compressibility 

(7)  Yield  Criteria 

(8)  Yield  zone  stress-strain  relations 

The  boundary  conditions  are  fixed  by  the  geometry  and  the  physical 
environment.   It  was  assumed  that  at  the  inner  bore  the  radial  stress 
was  zero  and  that  the  outer  surface  had  an  induced  pressure  caused  by 
the  restraining  effect  that  the  case  applied  to  the  propellant  as  the 
motor  was  cooled.   Plane  strain  was  the  assumed  end  condition. 

For  mathematical  simplicity  in  many  plasticity  solutions,  the  material 
has  been  assumed  to  have  the  same  compressibility  characteristics  in  the 
yielded  region  as  in  the  elastic  region.   Furthermore,  the  material 
has  been  taken  to  be  incompressible,  which  for  metals  can  introduce 
considerable  error.   However,  this  assumption  is  quite  valid  when  dealing 
with  solid  propellant,  as  Wood's  experimental  evidence  has  shown.   Thus, 
it  was  assumed  that  the  material  was  incompressible  in  both  the  dewetted 
and  the  undewetted  regions.   This  gave  considerable  mathematical 
simplification  in  that  Poisson's  ratio, jJ,   was  assumed  to  be  a  constant 
and  equal  to  one  half  and  that  the  sum  of  the  dewetted  strain  increments 
was  equal  to  zero. 


dff  +  dt'Q      +  d£^  =  0   ■  (7) 

This  assumption  permitted  the  simplification  of  the  yield  criterion 

and  the  stress-strain  laws  to  useable  forms. 

Two  yield  criteria  were  frequently  encountered.   The  von  Mises, 

where 

J.  =  -2-°  (8) 
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and  the  Tresca  yield  condition,  in  which 

The  Tresca  condition  offered  great  mathematical  simplicity  for  the 
cylinder  problem  facilitating  solutions  that  were  otherwise  out  of  reach, 
provided  that  the  relative  magnitudes  of  the  three  principle  stresses 
were  known  a  priori.   It  was  shown  that  in  this  problem,  the  axial  stress 
was  the  middle  stress,  thus  there  was  no  need  to  calculate  it.   Addi- 
tionally, its  contribution  to  plastic  yielding  was  zero. 

On  the  other  hand,  the  von  Mises  criterion  required  no  prior 
knowledge  of  the  stresses,  but  it  yielded  nonlinear  differential 
equations.   Due  to  the  fact  that  investigators  have  demonstrated 
experimentally  that  the  von  Mises  criterion  provides  a  closer  correla- 
tion with  a  description  of  yield  surface  than  the  Tresca  condition, 
solutions  have  been  sought  and  obtained  for  the  pressure  problem  using 
the  von  Mises  criterion.   The  thermal  problem,  however,  was  more  difficult 
in  that  a  thermal  term  was  added  to  the  stress-strain  equation.   The 
addition  prevented  a  simple  closed  form  solution  using  the  von  Mises 
criterion;  however,  it  had  no  effect  on  the  Tresca  condition.   Consequently, 
the  Tresca  condition,  even  though  less  accurate,  was  frequently  chosen 
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over  the  von  Mises  criterion.   Recall,  though,  that  from  dewetting 
theory  Wood  found  the  dewetting  criterion  to  be  of  the  von  Mises  form. 
It  was  therefore  desired  to  obtain  a  solution  using  the  von  Mises 
surface. 

There  are  two  plastic  flow  rules  that  have  been  used  exclusively  the 
Prandtl  Reuss  and  the  Hencky: 

Prandtl-Reuss 

where  d£„  is  the  equivalent  plastic  strain  and  oc  the  equivalent  stress. 
Hencky 

*  =  £  1<V  t <***•>} 

The  Prandtl-Reuss  equations  deal  with  incremental  theory  where  an 
increment  of  plastic  strain  is  so  small  that  it  can  be  considered  linear, 
The  total  plastic  strain  is  achieved  by  integration.   An  expansion  of 
equation  (5)  for  dewetting  materials  shows  it  to  be  identical  to  the 
Prandtl-Reuss  equations. 

Hencky  postulated  that  the  total  plastic  strain  would  have  the  same 
type  of  relationship  as  the  incremental  such  that  the  present  state  of 
plastic  strain  was  independent  of  its  loading  path.   This  is  not  true 
in  general,  and  the  Hencky  relationships  have  fallen  into  disuse. 
However,  under  special  circumstances  of  linear  strain  hardening  and  a 
monotomically  increasing  load,  the  two  flow  rules  coincide  (5). 
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Steele  presented  a  summary  of  the  combinations  of  these  parameters  that 
have  been  tried  by  different  researchers.   It  is  repeated  here  in  Table  I. 
Note  that  most  of  the  solutions  required  numerical  analysis.   Any 
numerical  solution  necessitates  a  loss  of  generality  and  makes  a  para- 
meter study  more  costly  and  difficult.   Thus  a  closed  form  solution  was 
desired.   By  making  simplifying  approximations  which  were  the  result 
of  either  the  special  geometry  assumed  or  the  unique  mechanical  properties 
of  solid  propellant,  it  was  possible  to  obtain  a  closed  form  solution 
to  the  thermal  problem  using  a  von  Mises  dewetting  surface. 
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III.   THEORETICAL  ANALYSIS 

Within  the  above  mentioned  framework  and  background,  the  theoretical 
analysis  was  undertaken.   Figure  2  describes  the  geometry  of  the  problem, 
and  the  governing  equations  are  given  below. 

Equilibrium 

dd'r   ^r-<£$    „  (12) 


dr 


=  0 


Strain  Compatability 

d^e    €g~  £f   =  0  (13) 

dr       r 

Stress-Strain 

E^ii^'^^V^*^  (14) 

where  £roj  are  the  dewetted  strain  components. 

The  problem  was  well  structured  in  that  there  were  sufficient 
equations  to  accommodate  the  unknowns;  thus,  a  unique  solution  could  be 
derived.   By  applying  the  plane  strain  assumption  and  the  incompressi- 
bility  condition  of  the  dewetted  zone,  the  last  of  equations  (14)  was 
solved  for  the  axial  stress 

6t   =  l/2(<Df+^e)-EotT  +•  E(  £^  +  el)  (15) 

Using  this  definition  of  the  axial  stress,  the  first  two  stress-strain 
relations  (equations  (14)  were  substituted  into  the  strain  compatability 
equation  (equation  (13))  yielding  a  differential  equation  in  terms  of 
oe  at*A  <£p .   The  tangential  stress,^,  was  eliminated  through  the  use 
of  the  equilibrium  equation  (equation  (12)),  yielding 
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-T 

(18) 


1  d^r  __  2   (£r-4)=  o 
E  dr     3      r 

Integrating  this  equation  twice  yielded  an  expression  for  the  radial 
stress. 

where  C,  and  C~  were  constants  of  integration.  Analytic  integration  is 
only  possible  when  Poisson's  ratio  is  one-half,  which  is  the  actual 
value  for  propellant. 

The  temperature  distribution  was  assumed  to  be  constant  as  previously 
mentioned.   Consequently,  the  temperature  integral  of  equation  (17)  was 
readily  evaluated. 

The  two  constants  of  integration  were  evaluated  by  applying  the  boundary 

conditions : 

at  r  =  a      <3r  =  0 

(19) 
at  r  =  b       <5r  =  p' 

where  p'  was  the  pressure  at  the  boundary  between  the 

propellant  and  the  case  and  was  caused  by  the  difference  between  the 

case's  coefficient  of  thermal  expansion  and  that. of  the  propellant. 

See  Appendix  A.   Application  of  the  first  boundary  condition  yielded 

C2  =  EC,  (20) 

Applying  the  second  boundary  condition  at  r  =  b,  gives 

Vjj&N-fJW*'*^)]  C21) 

Substituting  (19)  and  (21)  back  into  the  expression  for  6f  yielded 


(22) 


3  )*.  r 
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The  need  was  now  to  solve  the  remaining  dewetted  strain  integrals. 
As  mentioned  above,  it  had  been  assumed  that  the  Lindsey  dewetted 
stress  strain  equations  applied  and  that  solid  propellant  possessed 
linear  strain  hardening  characteristics.   Expanding  equations  (5)  gives 

As  previously  stated,  for  monotonically  increasing  loads  and  a  linearly 
strain  hardening  material,  the  Prandtl-Reuss  equations  degenerated  to 
the  Hencky  equations.   Within  this  approximation,  therefore,  it  was 
assumed  that  the  Lindsey  equations  took  the  form  of 

From  those  equations  it  is  seen  that 

fcr     fce      Z     £e  ^  «•     eJ  (25) 

=  ^[c£e-£^?vte1-£;tv(6;-£e^l',l 


(24) 


where 


£p 


6  =^  ~-p  [(<%^e^rV+^^] 


(26) 


The  axial  stress,    ,  was  replaced  by  equation  (15),  and  after  expanding 
the  terms  and  regrouping 


4=^[(<W-(S^-Tf 


(27) 

An  examination  for  order  of  magnitude  of  each  term  of  the  right  side 
of  equation  (27)  demonstrated  that  the  square  of  the  stress  difference 
clearly  dominated  the  entire  expression.   The  radial  and  tangential 
stresses  were  of  opposite  sign  and  such  size  that  their  difference  was 
relatively  large.   The  low  modulus  of  the  propellant  relative  to  that 
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of  some  other  material  such  as  steel  (see  Figure  13)  made  the  E^Tterm 
small  compared  to  the  stress  difference.   Additionally,  from  the  elastic 
solution,  it  is  seen  that  the  radial  strain  is  negative  while  the 
tangential  strain  is  positive.  Moreover,  the  absolute  magnitude  of  the 
radial  strain  is  larger  than  that  of  the  tangential  strain;  consequently 
the  sum  of  the  tangential  and  radial  dewetted  strain  components  must  be 
negative,  assuming  that  the  direction  of  the  dewetted  strain  increments 
are  in  the  same  direction  as  their  elastic  components.   So  it  seems 
that  the  already  small  EdX  term  is  reduced  even  more  by  the  addition  of 
this  negative  quantity.   Thus,  the  second  term  inside  the  radial  was 
neglected.   The  expression  for  the  equivalent  stress  became 

<£«  =J3K-«el  (28) 

which  simplification  made  possible  a  closed  form  solution  without 
appreciable  loss  in  accuracy.   The  terms  neglected  were  estimated  to 
constitute  at  most  a  2%  error. 

Using  equation  (28) ,  equation  (25)  reduced  to 

t'x-  *e  =  \  Z?  S3n(V<Q)  (29) 

where  Sffft (<V ^^  <  O    for  the  thermal  problem. 

The  significance  of  assuming  linear  strain  hardening   actually 
meant   that  the  stress-strain  curve  was  approximated  by  two  straight 
lines,  one  having  the  slope  of  E  and  the  second  having  the  slope  of 
mE  as  shown  in  Figure  4.   From  the  figure  an  expression  for  the  dewetted 
region  was  found  to  be 

By  taking  equation  (28)  and  applying  the  equilibrium  equation  (12)  it 
is  seen  that 

<^JT  ^-r^  (31) 
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Using  equations  (31,  30,  29)  permitted  direct  integration  of  the 
dewetted  strain  integral  in  equation  (22)  such  that  the  radial  stress 
expression  became 

A  =  -m_  LorV)^  i-   <_  ^iiL^r^-4^0U(a/r^]       02) 

Two  new  unknown  constants  appeared  in  this  equation,  Or  and  r  . 
The  critical  radius,  r  ,  was  the  radius  of  the  boundary  between  the 
dewetted  and  the  elastic  zones,  and  or  was  the  radial  stress  at  this 
boundary.   Both  of  these  were  constant  and  were  functions  of  the  thermal 
loading  only.  Gr  was  related  to  r  by  setting  r  equal  to  r  in  equation 
(32).   After  substituting  for  p' 

4=   ftf-a^t  aft     II     1_J~±  *  (33) 

where  p'  and  D  came  from  the  elastic  solution.   (See  Appendix  A.) 

It  was  noted  that  at  the  boundary  between  the  dewetted  and  undewetted 

zones,  the  radial  stress  must  be  continuous  across  the  boundary. 

Further,  note  at  that  location  that  the  dewetted  condition  must  also  be 

met  by  the  stresses  in  both  the  elastic  and  the  dewetted  regions. 
The  yield  condition,  equation  (1)  became 

(6r-4e-)    =^B0V<^ff  (34) 

The  second  term  under  the  radical  was  considered  negligible  for  reasons 

already  given.   The  expression  for  or  and  <Se  were  taken  from  Appendix  A. 

The  combined  loading  case  was  used,  and  the  radius  was  set  equal  to  the 


critical  radius  (r=r  ) ,  such  that 
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Combining  equations  (34)  and  (35)  and  solving  for  r  yielded 


t  = 


3  bliii*^       tf(m'V^,A  -r 


z 


(36) 


6tK 

It  would  have  been  possible  to  solve  equation  (34)  and  (35)  for  O-.  in 

terms  of  an  unknown  rc  and  to  set  this  new  expression  equal  to  equation 
(33),  solving  for  r  directly.   However,  such  a  course  of  action  would 
have  resulted  in  a  transcendental,  tenth  order  expression  and  would  have 
been  of  little  practical  value.   A  simpler  method  was  to  iterate  with 
equations  (36  and  (33) .   Assume  a  value  for  6,  (say  or  =  0)  and  solve  for 
r  by  using  equation  (36).   This  value  for  r  was  used  in  equation  (33) 
to  solve  for  a  better  value  of  Of  •   This  process  was  repeated  until 
the  resulting  values  of  r  converged. 

All  unknown  constants  of  equation  (32)  were  known  and  the  complete 
expression  for  the  radial  stress  was  obtained: 

E1S  6r=^V^{-f'-^v^V^l^^w=]     o7) 


4=-^^'-^^^}^^^ 


(38) 
Through  the  application  of  the  equilibrium  equation,  the  tangential 


stress  was  obtained  directly. 


(39) 


The  strains  were  obtained  by  application  of  the  stress-strain 
equations. 

r    H        m    E     ?-  HE>w\  (A0) 

3  _,T   .   <^  <£> 
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IV.   EXPERIMENTAL  RESULTS 

There  was  a  dual  goal  in  the  experimental  work  done.   First,  it 
was  necessary  to  experimentally  determine  the  constant,  B(T)  in  equation 
(2),  for  the  particular  PBAN  propellant  used.   This  was  accomplished  by 
conducting  uniaxial  tensile  tests  over  a  range  of  temperatures  in  the 
gas  dilatometer  of  United  Technology  Center,  Sunnyvale,  California. 
For  uniaxial  tension,  Equation  (2)  reduced  to: 

B(T)  =   ^  (41) 

The  dilatometer  (see  Figures  5,  6,  7,  8)  utilized  the  Perfect  Gas 
Law  as  its  basic  operational  principle.   The  specimen  was  mounted  within 
a  sealed  chamber  and  pull-rods  extended  from  either  end  of  the  sample 
through  the  ends  of  the  chamber  and  attached  to  the  cross-head  of  an 
Instron  testing  machine.   Any  change  in  the  volume  of  the  sample  while 
under  load  resulted  in  a  change  in  the  chamber's  pressure  creating  a 
pressure  difference  across  a  diaphram  separating  the  test  chamber  from 
an  adjoining  reference  cavity.   A  baratron  pressure  transducer  converted 
this  pressure  difference  into  an  electrical  signal,  which  was  fed  to 
one  pen  of  a  two-pen  X-Y  plotter.   The  signal  for  the  second  pen  came 
from  a  load  cell  in  the  cross-head  of  the  Instron  machine.   The  raw  data 
thus  obtained  were  in  the  form  of  two  superimposed  plots;  one  was  the  time 
history  of  the  load,  and  the  second  was  the  time  history  of  the  dilata- 
tion.  By  using  the  constant  cross-head  speed  and  premeasuring  both 
the  cross  sectional  area  and  gage  length,  the  data  was  reduced  to  plots 
of  stress  vs  strain  and  dilatation  vs  strain.   (See  Figure  9) 
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Changes  in  the  ambient  temperature  caused  a  problem  since  it 
affected  the  readings  of  the  pressure  transducer  through  the  same  Perfect 
Gas  Law.   In  order  to  isolate  the  experiment  from  ambient  temperature 
fluctuations  and  also  provide  a  method  of  controlling  the  temperature 
at  which  the  experiment  was  conducted,  the  dilatometer  was  enclosed  in 
a  temperature  box,  which  possessed  the  capability  of  maintaining  a 
preselected  temperature  to  within  a  degree  over  the  duration  of  the 
experiment.   Although  the  ambient  temperature  within  the  temperature 
box  was  fairly  constant,  at  subzero  temperatures,  heat  leaked  into  the 
dilatometer  probably  by  conduction  through  the  push  rods.   A  large 
temperature  drift  would  develop  after  about  one  half  hour  and  this 
dictated  the  maximum  duration  of  any  individual  run. 

The  second  goal  was  to  experimentally  determine  the  mechanical 
properties  of  the  subject  propellant  as  functions  of  temperature; 
namely,  the  uniaxial  dewetting  stress,  o   »  and  the  elastic  modules,  E 
and  m.   (See  Figures  10,  11,  12) 

The  test  samples  were  square  with  sides  machined  to  approximately 
one  half  inch,  and  they  had  a  gage  length  of  approximately  2  3/4  inches. 
Redwood  tabs  were  bonded  to  each  end  of  the  sample  to  make  hardware 
attachment  possible.   Tests  were  conducted  at  21°C,  3°C,  -9°C,  -21°C, 
and  -52°C,  with  four  tests  being  conducted  at  the  three  higher  tempera- 
tures, two  tests  conducted  at  -21°C  and  one  test  conducted  at  -52°C. 
Experimental  difficulties  with  the  equipment  at  the  subzero  temperatures 
resulted  in  the  reduced  number  of  runs  at  the  lower  temperature. 
Experimental  scatter  of  seven  to  ten  percent  was  found  in  the  dilatation 
and  stress-strain  curves.   However,  only  about  five  percent  scatter  was 
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found  among  the  parameters  measured  (i.e.,  E,  m,  and  oo  ) .   Due  to  the 
fragile  nature  of  solid  propellant  specimens,  this  was  considered  quite 
reasonable  results. 

After  reducing  the  raw  data  to  the  stress  vs  strain  and  dilatation 
vs  strain  plots  mentioned  above,  the  problem  of  identifying  6  developed. 
Three  methods  were  available  from  the  literature:   (1)  interpreting  the 
stress  which  corresponds  to  the  lift-off  point  of  the  dilatation  curve 
as  od  .   This  method  provided  results  which  were  quite  repeatable; 
however,  the  stress  strain  curve  did  not  depart  from  its  initial  linear 
relationship  at  this  point.   (2)  The  dilatation  curves  eventually 
assume  a  constant  slope.   Physically  the  total  amount  of  solid  material 
which  could  dewet  had  dewetted  at  this  point  and  further  dilatation 
resulted  only  from  increasing  the  size  of  previously  formed  voids.   The 
stress  corresponding  to  the  extrapolation  of  this  constant  slope  line 
back  to  the  zero  stress  point  (pt.  B), of  Figure  13  would  be  considred  as 
the  yield  stress,  o  .   This  could  be  physically  interpreted  as  assuming 

0 

that  no  dewetting  occurred  until  the  yield  stress  was  reached;  then  all 
the  dewetting  occurred  simultaneously.   This  was  tantamount  to  assuming 
that  the  propellant  was  perfectly  plastic.   However,  at  the  lower 
temperatures,  specimen  failure  resulted  before  such  a  linear  portion  to 
the  dilatation  curve  was  reached.   (3)  The  third  method  was  to  take  the 
stress  at  the  point  where  the  stress-strain  curve  initially  departed 
from  a  linear  relationship.   This  method  reflected  the  linear-strain 
hardening  assumption,  and  provided  consistent  reproduceable  results 
from  sample  to  sample.   It  neglected  the  dilatation  curve  altogether. 
This  last  method  was  adopted  primarily  because  it  met  the  condition 
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of  the  linear  strain  hardening  assumption.  A  comparison  of  these 
three  methods  was  given  in  Table  II,  where  it  is  seen  that  methods 
(2)  and  (3)  yield  similar  results. 

The  test  speeds  were  chosen  as  slow  as  possible,  yet  fracture  had  to 
be  reached  within  the  time  limit  mentioned  above.   At  the  lower  tempera- 
tures, it  was  assumed  that  the  speeds  used  maintained  the  specimens  in 
equilibrium  and  that  the  load  did  not  produce  viscoelastic  effects. 

In  summary,  the  yield  stress,  <j  ,  was  taken  as  the  stress  at  which 
the  average  stress  strain  curve  first  departed  from  linearity  for  all 
runs  at  any  one  temperature.   The  modulus,  E,  was  the  average  of  the 
elastic  slopes  of  the  stress-strain  diagrams  for  any  given  temperature. 
The  strain  hardening  exponent,  m,  was  the  ratio  of  the  slope  for  the 
least  squares  straight  line  fit  of  the  non-linear  portion  of  the  stress 
strain  curve,  to  the  initial  modulus. 
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V.   RESULTS  AND  CONCLUSIONS 

A.   RESULTS 

In  order  to  compare  the  results  predicted  by  dewetting  theory  with 

those  obtained  from  linear  elastic  theory,  a  sample  problem  was  worked. 

A  PBAN  propellant  of  current  interest  was  used,  and  its  mechanical 

properties  were  experimentally  determined,  (Figures  11,  12).   The  other 

required  parameters  were  taken  from  Appendix  II  of  Reference  (4) . 

For  the  propellant:    V=   1/2  c*  =  2.21xl0~5  in./in./°F 
For  the  steel  case:    1^  =  1/3  dc=  .65xl0~5   in./inc./°F 

E  =  3xl07psi 
c        r 

Geometrical  parameters:  a=  3in.  b=  12in.  h=0.015in. 

The  rocket  motor  was  cooled  from  160CF  to  -60°F,  and  both  elastic 

and  dewetted  stress-strain  analyses  were  made.   Linear  strain  hardening 

was  assumed  which  enabled  calculations  to  be  performed  in  two  steps: 

one  elastic  step  and  one  dewetted  step.   The  temperature  was  not 

incrementally  stepped  down.   Equations  (32)  and  (36)  were  iterated  in 

order  to  find  the  critical  radius,  r  ,  and  the  interface  stress, <Sr  . 

c'  c 

These  results  were  plotted  in  Figure  (14).   It  was  interesting  to  note 
that  the  growth  of  the  yield  surface  was  linear  with  temperature. 
Furthermore,  dewetting  did  not  begin  until  the  temperature  had  been 
decreased  to  10°F,  and  at  the  lower  temperature  limit  the  yield  surface 
had  penetrated  to  a  radius  of  only  9  1/2  inches.   The  interface  pressure 
also  increased  linearly  as  the  temperature  was  lowered  reflecting  the 
influence  of  temperature  on  both  the  yield  criterion  and  the  propellant' s 
mechanical  properties. 
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The  greatest  effect  of  dewetting  was  upon  the  tangential  and  radial 
stresses.   (Figures  15,  16)   The  tangential  stress  at  the  bore,  for  both 
the  elastic  and  the  dewetting  theories,  was  plotted  against  temperature. 
It  was  seen  that  dewetting  substantially  relieved  the  stress.   The 
radial  stress  was  plotted  across  the  motor's  thickness  for  a  characteris- 
tic temperature  of  -8.0°F.   Again  it  was  shown  that  stress  was  relieved 
within  the  dewetted  region.   The  radial  stress  then  approached  the 
elastic  solution  as  one  moved  deeper  into  the  undewetted  elastic 
material  and  away  from  the  critical  radius.   However,  a  pronounced 
effect  of  the  induced  pressure  at  the  inner  radius  of  the  undewetted 
region  was  still  evident. 

Initially,  the  behavior  of  the  tangential  strain  at  the  bore  was 
quite  surprising.   (See  Figure  (17))  It  had  been  anticipated  that  this 
strain  would  have  been  negative  and  would  have  increased  in  magnitude 
as  the  temperature  was  lowered.   However  upon  looking  at  the  elastic 
tangential  strain  equation  (A-l) ,  one  sees  that  it  consists  of  two  terms. 
The  first  term  represents  the  mechanical  strain  imposed  by  the  case 
upon  the  propellant.   It  is  positive  and  dominates  the  second  term, 
the  thermal  strain,  which  is  negative  in  the  cool-down  problem.   However, 
as  the  temperature  is  lowered  the  mechanical  properties  of  the  pro- 
pellant change  significantly,  and  the  negative  T  term  grows  in  relative 
magnitude,  actually  forcing  the  total  strain  to  decrease.   Thus,  it  was 
observed  that  the  maximum  strain  occurred  at  about  A0°F  instead  of  at 
the  lower  temperature  extremum. 
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B.   CONCLUSIONS 

A  point  worthy  of  note  is  the  curve  fit  of  the  m  curve  coefficient 
of  linear  strain  hardening  in  figure  12.   Equipment  limitations  at  the 
subzero  temperatures,  which  restricted  multiple  runs  decreased  confidence 
in  the  values  of  m.   Consequently,  the  specific  mechanical  properties 
measured  at  -5.8°F  and  -60°F  were  used  to  calculate  the  bore  hoop 
strain  at  those  two  temperatures,  and  they  appear  on  the  plot  in  Figure  14 
which  uses  the  curve  fit  for  m.   It  was  concluded  that  a  curve  fit  which 
averaged  these  values  provided,  a  smoother  extension  of  the  elastic 
solution  than  a  curve  which  connected  each  data  point.   Further  tests 
would  prove  useful  to  more  definitely  fix  the  m  curve. 

From  a  comparison  of  the  elastic  and  the  dewetting  solutions  for 
both  stresses  and  strains  it  was  concluded  that  dewetting  has  little 
effect  on  the  strains.   In  an  incompressible  material  the  strains  are 
fixed  by  geometrical  constraints  and  not  by  mechanical  properties. 
However,  stresses  are  tied  to  the  mechanical  properties  and  thus  are 
greatly  affected  by  dewetting. 

« 

It  was  further  concluded  that  a  thermo-dewetting  solution  could 
be  obtained  in  a  relatively  simple  form  and  that  better  results  for  stress 
would  be  obtained  by  using  dewetting  theory  vice  elastic  theory. 
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TABLE  II 
TABULATION  OF  EXPERIMENTAL  DATA 


Temp 

cr 

o 

E 

m 

21° 

94 

806 

,417 

3.5 

118 

983 

.545 

-9 

137 

1600 

.49 

-2L 

159 

8840 

.369 

-52,5 

231 

23100 

.539 
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TABLE  III, 


COMPARISON  OF  METHODS  OF  READING 
TENSILE  YIELD  STRESS 


TEMPERATURE 
(°C) 


22° 


3.5 

•9 
■21 
■  52.5 


I 

<£  (psi) 
II 

III 

46 

84.5 

94 

72.5 

107 

118 

105 

• 

137 

137 

• 

159 

201 

* 

231 

^Dilatation  Curve  did  not  reach  constant  slope. 

I„   o^  was  stress  corresponding  to  the  lift  off  point 
of  the  Dilatation  Curve. 

II.    69   was  stress  corresponding  to  extrapolation  of 

of  the  constant  Dilatation  Curve  back  to  the  strain 
axis  o 

III.    60   was  stress  at  which  stress-strain  curve  initially 
departed  from  a  linear  relationship. 


28 


Dewetted  Propellant 


Rocket  Motor  Geometry 
Figure  1  . 
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Circular 
Cross- section 


Dewetted  Yield  Surface  in  Principle  Stress  Space 

Figure  2.  . 
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Intersection  of  Thermal  Stress  Load 
with  Dewetting  Surface 

Figure  3. 
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Linear  Strain  Hardening 
Figure  4. 
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Figure    5 
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Figure    7 
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APPENDIX  A 

ELASTIC  SOLUTION 

Williams,  Blatz  and  Schapery  (7)  of  the  Guggenheim  Aeronautical 
Laboratory  (California  Institute  of  Technology)  analyzed  solid  pro- 
pellant  rocket  motors  of  the  same  geometry  under  many  different  loading 
modes.   Their  basic  assumptions  were  that 

1.  the  propellant  had  material  properties  which  were 

a.  Isotropic 

b.  Homogeneous 
d.   Continuous 

2.  the  strains  were  small. 

Their  solution  to  the  thermally  loaded,  cylindrical,  rocket  motor, 
assuming  a  condition  of  plane  strain,  was: 

A  .    \>Vp«o-0H  «»frg>  ft.a,  -  r%JL* 


(v-a.'i        (yiNff*f\  V       (        k 


(A-l) 


where 





46 


And  for  the  case  of  an  internally  pressurized  rocket  motor  their 
solution  was : 

r"     r  Kb*- oft         v>l-a> 
G*       rKkl-a?-  tf-fl> 

3  a?  ,fx 
'         D 


where  p.  =  internal  pressure 
1 


In  considering  the  effects  of  dewetting,  the  inner  boundary  condition 
changed.   The  inner  radius  of  the  elastic  region  was  the  same  as  the 
critical  radius,  i.e.,  it  was  equal  to  the  radius  of  the  dewetting 
surface.   As  the  dewetting  surface  moved  outward,  the  inner  radius 
of  the  elastic  zone  also  increased;  furthermore  the  inner  boundary  was 
no  longer  stress  free.   From  continuity  of  the  radial  stress  across  the 
dewetting  surface,  it  was  seen  that  the  inner  boundary  now  had  a  stress 
equal  toor  induced  upon  it.   This  elastic  solution  was  obtained  by  com- 
bining A-l  and  A-2,  and  setting  p^  equal  to <£r  ,  and  the  inner  radius, 
a,  equal  to  the  critical  radius,  rc 


4= 


<Sr,  *~c  1"°-  ;(tf-V^      *T>        ___     (A-3) 
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